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The calculation of the NNLO QCD corrections to the partonic process gg — > H is 
outlined. For the coupling of the Higgs boson to the gluons we use an effective Lagrangian 
in the limit of a heavy top quark. The focus is on the evaluation of the virtual two-loop 
corrections. It is shown that the leading pole terms are in agreement with the general 
formula by Catani. 

1. Introduction 

Gluon fusion will be the dominant production mechanism of a Standard Model 
Higgs boson at the LHC. For a Higgs mass between 100 and 200 GeV, gluon fusion 
exceeds all other production channels by a factor ranging from five to eightLL QCD 
radiative corrections to this process are found to be more than 50% at NLOa. The 
evaluation of the NNLO corrections (0(a^)) is therefore highly desirable. 

The first step towards the full result was taken recently!, when the virtual NNLO 
corrections were evaluated. Although infra-red divergent, the result could be used 
to deduce an expectation of the size of the full NNLO corrections. This estimate 
turns out to be roughly 10-20%, indicating that the perturbative expansion is valid. 

There are several more steps to be taken in order to arrive at a prediction of 
the inclusive cross section at NNLO. First, single and double real radiation have 
to be added to the virtual corrections and all terms must be renormalized. After 
factorization of the soft singularities into the splitting functions, this leads to an 
IR and UV finite result. Finally, the latter has to be folded with the corresponding 
parton distribution functions at NNLO whose evaluation is still awaited. 

In the following we will focus mainly on the virtual corrections, in particular 
on the comparison of the leading poles in e = (4 — D)/2 with a general result by 
Catani! (D is the space-time dimension which is used for the regularization of the 
integrals). For more details on the actual calculation and the results we refer to 
Refl 

2. Virtual two-loop corrections 

It has been shown! that the limit of a heavy top quark is well justified in the 
process gg — > H . Technically this means that the top quark loop that mediates the 
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2 Higgs production in gluon fusion to O(a^) 



coupling of the gluons to the Higgs boson can be replaced by an effective vertex. 
In this way, the NNLO contribution is represented by two-loop vertex diagrams 
with two massless on-shell legs. Such diagrams can be evaluated by mapping them 
onto three-loop massless propagator diagramsu and reducing them through the well- 
known integration-by-parts algorithm!]. 

We write the virtual contribution to the cross section for the process gg — > H as a 

. „. 9 2 



^virt — 



47T 

'} 1 — 



Cx(a s ) 



5(1 -z) 



l-(3{a s )/e) 256(1 - e) 



2L(D 



a {2) + 0(a 3 s ) 



v 2 [ ' 256(1 -e) 



1 + 

7T 



/?(a.) = - (^i) 2 /3l+ 0( a 3), ft 

7T \ TT / 



(£) «& + 0(a») 
33 - 2n z 



12 



153 - 19n; 
24 



(1) 



Ci(a s ) is the coefficient function for the <?gif vertex and contains the residual 
logarithmic top mass dependence!] (l t = ln(/i, 2 /M 2 ), with M t the on-shell top quark 
mass, ni = 5 denotes the number of light quark flavors, and a s = ai°\{j, 2 ) is the 
running coupling for five active flavors) 
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According to Catani's general resultu, the UV-renormalized NNLO contribution 
can be cast into the following form: 



aW=ll W (e)a£l+ 1 -lW(e 
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where aj^ is finite as e 
J (1) (e) reads 6 ( 7E = 0.577216 . . . ) 



0. In our case of two incoming gluons, the operator 
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"Note that Eq. (9) of Ref|f erroneously has an additional factor of Mj : . 
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Note that due to a misprint (A^ = +1 instead of \ 



-1) the original formulaQ results in 



a different sign of the unitary phase. A clarifying conversation on this issue with S. Catani and 
S. Dittmaier is acknowledged. 
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H^ 2 \e) contains only single poles in e. Thus, Eqs. (0)-(|5|) determine the poles of 
order 1/e 4 , 1/e 3 , and 1/e 2 of the NNLO amplitude a^. They can now be compared 
with the explicit evaluation of the Feynman diagrams at NNLOu, and full agreement 
is found? This provides a non-trivial check on the results of Reffl. From the latter, 
one may now extract the undetermined pieces (e) of Eq. (||) and of Eq. (^). 
Their values depend on whether or not one keeps the higher order terms in e when 
expanding the exponentials and the T-functions in Eqs. (d)-(H). Thus, in order to 
prevent confusion, we refrain from giving explicit expressions for H^ 2 \e) and a|^. 
Instead, we advise the interested reader to use the results of Ref.i. 

3. Real radiation 

Let us give a brief description of the contributions coming from the real emission 
of quarks and gluons. The divergences of the virtual two-loop corrections will be 
canceled by the soft contributions of the following processes: gg — > Hg to one-loop 
order, and gg — > Hgg, gg —* Hqq at tree level (assuming the effective ggH vertex). 
After adding these contributions to the virtual corrections, there are still infra-red 
singularities left. These have to be absorbed by the mass factorization procedure, 
using the Altarelli-Parisi splitting functions at NNLO. 

Explicit results for the full partonic calculation will be given elsewhereu. For 
the contribution oc 5(1 — s/Mg) (s is the cms energy of the gg-system) we find a 
NNLO correction of order 10%, in agreement with the numerical estimate of Refit! 
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c It is important to note that the prescriptions for UV-renormalization of RefQ and Ref£| differ in 
the treatment of terms oc 7e and oc ln47r. For convenience, we apply the former prescription here. 
For infra-red finite quantities, the two prescriptions are equivalent up to terms of order e. 



